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Lattice Model of Sweeping Interface for Drying Process in Water-Granule Mixture 
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Based on the invasion percolation model, a lattice model for the sweeping interface dynamics 
is constructed to describe the pattern forming process by a sweeping interface upon drying the 
water-granule mixture. The model is shown to produce labyrinthine patterns similar to those 
found in the experiment [Yamazaki and Mizuguchi, J. Phys. Soc. Jpn. 69 (2000) 2387]. Upon 
changing the initial granular density, resulting patterns undergo the percolation transition, but 
estimated critical exponents are different from those of the conventional percolation. Loopless 
structure of clusters in the patterns produced by the sweeping dynamics seems to influence the 
nature of the transition. 

KEYWORDS: interface dynamics, sweeping interface, pattern formation, invasion percolation, granular 
system, drying process 



1. Introduction 

During drying process, it often happens that mate- 
rial dispersed in fluid condenses and leaves patterns on 
a substrate. 1 A typical one is dry stain made by cof- 
fee spilt over a table. Yamazaki and Mizuguchi demon- 
strated 2 that the water-granule mixture that is confined 
in the narrow gap between two glass plates produces 
much more complicated patterns than coffee stain upon 
drying; Labyrinthine patterns of granules emerge as gran- 
ules are swept by one-dimensional water-air interfaces as 
they recede during the drying process. 

Physical processes are as follows; As the water evapo- 
rates from the gap between the plates, the volume of the 
wet region tends to shrink and the pressure in the wa- 
ter decreases, which gives the driving force to move the 
interface between the water and the air; Although the 
pressure difference should be almost uniform over the 
whole interface, the interface does not move uniformly, 
but only the weakest part moves at a time. As the in- 
terface moves, it sweeps the dispersed granules to collect 
them along it; Consequently, the interface becomes more 
difficult to move due to the friction of granules with the 
glass plates. At some parts of the interface, the granular 
density eventually exceeds the threshold where granules 
get stuck to form a pattern after drying. 

For this phenomenon, several models have been pro- 
posed. In their original paper on the experiment, Ya- 
mazaki and Mizuguchi analyzed the elementary process 
to estimate the thickness of the granule region assuming 
the granular region as elastic material. 2 A phase field 
model with the granular density field 3,4 and the bound- 
ary dynamics 5 have been constructed, based on the ob- 
servation that the sweeping phenomenon can be regarded 
as the small diffusion limit of the crystal growth. 6 ' 7 

In this paper, focusing on similarities in dynamics be- 
tween the sweeping process described above and the in- 
vasion percolation, 8 we present a simple lattice model of 
sweeping interface for pattern formation. We performed 
numerical simulation on the model and analyzed the re- 
sulting patterns. 



2. Lattice Model of Sweeping Interface 

The model consists of two variables at each site of the 
triangular lattice: the state variable Si and the granular 
density fi on the site i. The state variable Si takes the 
values either or 1 depending upon whether the site i is 
dry or wet: 



(if the site i is dry) 

1 (if the site i is wet) 



(1) 



The granular density fi takes zero or a positive real num- 
ber representing the quantity of granules in the i'th cell. 

Initially, all the sites in the system are wet and the 
granules are distributed uniformly with some random- 
ness, thus the initial value of the state variables Sj's are 
1 for all the sites, and the granule density fa's are given 
random values from the uniform distribution with the 
width A/ and the upper bound f M : € [f M - A/, f M }. 

The interface site is defined as a wet site adjacent to 
a dry site. We assume the sites outside of the system 
are dry, thus the wet regions are surrounded by the one 
dimensional chains of the interface sites, which we call 
the interface. 

Upon drying, the water volume shrinks, then the pres- 
sure decreases in the water region. The pressure differ- 
ence between the wet and the dry regions pushes the 
interface toward the wet region, but the site that actu- 
ally moves is the least resistive site. The resistance, or 
the strength of the interface sites, against the interface 
motion is determined by the granular density fi and the 
local configuration of interface; The granules causes fric- 
tion with the glass plates, while the interface sites sur- 
rounded by the dry sites are easier to dry due to the 
surface tension that tends to make the interface straight. 
In order to take these effects in a simple way, we intro- 
duce the strength r% of the interface site i as a function 
of its granular density fi and the number of neighboring 
dry sites n». We employ the simple form 



fi if rii < n s 

fi - 7 if i%i > n s 



(2) 
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with two parameters 7 and n Sl which represent the sur- 
face tension effect. 

Sweeping of granules takes place when the interface 
moves; Upon drying a wet interface site, the granules on 
this dried site are swept away to the neighboring wet 
sites, as long as the granular density of the site is not 
high enough for the granules to get stuck. If the granule 
density exceeds a threshold value g t h, the granules get 
stuck between the glass plates and cannot move, thus 
they stick to the sites and form a pattern after the whole 
system is dried. 

To represent the above processes, the dynamics for the 
lattice model of the sweeping interface is defined as fol- 
lows, (i) Out of the sites with ft < g t h, pick the weakest 
interface site i, whose strength is smallest. If all the 
interface sites have /j > 5th, then pick the weakest site 
from them, (ii) Dry the site i by changing its state vari- 
able Si from 1 to 0. (iii) Redistribute fi to the neighboring 
wet sites equally if fi is smaller than g t ^ and there exist 
some wet neighbor sites; Otherwise do nothing. Repeat 
the steps (i) - (iii) until all the sites becomes dry. 

There are five parameters that characterize the model: 
Af and /m to characterize the initial distribution of 
granular density, 7 and n s to control the effect of sur- 
face tension on the interface motion, and gth, that is the 
upper limit for the sweeping to take place. Out of four pa- 
rameters, Af, Jm , 7, and g t h, we can eliminate one that 
sets the scale for the granular density; We take g t h = 1- 

In the following simulations, we fix the parameter Af 
to be 0.5/m for simplicity, and take n s = 3, i.e. half of 
the number of neighboring sites in the triangular lattice 
that we are using. Then we have only two parameters, 
/m and 7. 

3. System behavior 

3.1 Patterns and surface tension effects 

First, we compare the patterns for 7/gth = with 
those for 7/gth = 1 in order to examine the surface ten- 
sion effects in Fig.l, where the sites with granules are 
represented by the solid circles whose sizes are propor- 
tional to the density. The initial grain densities are given 
by fivi/gth = 0.1, 0.5, and 0.9, and the system size Lx L 
with L = 64. One can see the surface tension effects es- 
pecially in the lower density cases; The grains tend to 
align in the systems with the surface tension j/gth = 1 
for /m = 0.1. These correspond to the patterns obtained 
in the experiment. 2 In the higher density cases, the sur- 
face tension effect is small. In the rest of the paper, we 
study only the case of 7/gth = 1 and change /m- 

Note that, in the case oi"f/g t h = 1, the pattern consists 
of winding paths delineated by walls of granular sites; 
the width of the paths are proportional to 1//m, while 
that of the walls is always of order of one lattice site 
because we assume that granules are swept only to the 
neighboring sites in the present model. 

3.2 Invading process 

Fig. 2 shows time development of the drying process 
for fM/gth = 0.871 for the system with L = 64. While 
one site is dried at each time step, the spatial develop- 
ment is not homogeneous but intermittent as in the case 




Fig. 1. Patterns produced by the sweeping interface model with 
and without the surface tension effect (7 = for the left column 
and -f/gth = 1 f° r the right column). The sites with the granules 
are represented by the solid circles whose sizes are proportional 
with the granule density. The initial granular density is given by 
/m = 0.1(the top row), 0.5(thc middle row), and 0.9(thc bottom 
row) with Af = 0.5/m and ri s = 3. The linear system size is 
L = 64. 

of the invasion percolation model. The difference in the 
pattern from the invasion percolation model is that the 
invaded (or dried) region develops a winding path struc- 
ture, and the basic mode of development is extending 
narrow paths. This is due to the sweeping process, by 
which the granules on a dried site are swept away to the 
adjacent wet sites; This makes walls on both sides of the 
path, consequently, it is easier to advance ahead making 
a narrow path than to widen the swept region beyond 
a certain width. This path extending process leads to a 
labyrinthine pattern of granules when the drying process 
is completed. 

Both of the two plates in the lower row of Fig. 2 rep- 
resent the same pattern obtained by this process; The 
lower left plate shows the grain density distribution by 
solid circles whose size is proportional to the density. 
The pattern consists of strings of the sites with grains. 
If we define a cluster as a set of connected sites with 
non-zero grain density, the whole pattern can be decom- 
posed into clusters; In the lower right plate of Fig. 2, its 
cluster structure is shown by coloring. The clusters are 
of highly branched structure, but there are virtually no 
loops of granular sites; The shape of clusters is compact 
and their fractal dimension appears to be very close to 
two. One can see that one of the clusters is percolating 
from the top to the bottom; In this case, we call the 
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Fig. 2. (Color online) Time development and final patterns. The 
upper two plates show the drying process. The blue regions de- 
note the wet regions. The lower left plate shows the final pattern 
after drying with the granule density represented by the solid 
circles. The lower right plate shows the cluster structure of the 
same pattern. The parameters are $m 1 9th — 0.871, A/ = 0.5/mi 
n s = 3, ■y/gth = 1 and L = 64. 



system percolated. 

4. Percolation transition of the sweeping inter- 
face model 

In the following, we will examine the percolation tran- 
sition in the pattern produced by the sweeping interface 
model. 

4-1 Occupation ratio 

The proportion p of the sites with grains after drying 
is not an input parameter but is obtained by perform- 
ing simulations with the initial distribution given by /m- 
In Fig. 3(a), p is plotted against fu/gth- The plot is for 
the system with L — 256, but the size dependence is 
small when L is larger. The ratio p is an almost linearly 
increasing function of /w/ffth) but one can see a small 
bump near /Af/ffth ~ 0.9, orp« 0.5; This roughly cor- 
responds with the percolation transition point. We will 
not pursue this any further in this paper although we do 
not understand its nature yet. 

4-2 Percolation probability 

Fig. 3(b) shows the percolation probability II as a func- 
tion of p for various lattice size L. The lines show fitting 
curves by the error function as 



n( P ;L)=Erf((p-p B0 (L))/Ap(L) 



where the error function is defined by 



Erf(x) = -= 

\/7T 



dt: 



(3) 



(4) 



P5o(L) is the occupied site proportion where II = 0.5, 
and Ap(L) is the width of the percolation transition for 
the finite system of the size L. For lower granular den- 
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Fig. 3. (a) Proportion of occupied sites p vs. Jm for the system 
size L = 256. (b) The percolation probability II vs. the occupa- 
tion proportion p for various system sizes L. The system sizes 
are L =124, 256, 512, 1024, 2048, 4096, and 8192 from the left 
curve to the right. The lines denote the fitting function by the 
error functions with pso(L) and Ap (see text). Each plot repre- 
sents average over 200 ~ 1000 realizations. The vertical line is 
drawn at p = p c = 0.5221. (c) pso(L) — Pc and Ap vs. L in the 
logarithmic scale. p^o(L) — p c is plotted with p c = 0.5221 in the 
left scale and Ap in the right scale. The solid lines represent the 
lines with the slope —0.50. The other parameters of the model 
are A/ = 0.5/ M , l/9th = b and n s = 3. 



sity, the system is not percolated, while the percolation 
probability is almost one at higher density; The curve is 
steeper, or Ap(L) is smaller, for the larger system, which 
suggests that there exists a sharp transition at a finite 
density in the infinite system size limit as in the case of 
the conventional percolation problem. 

In comparison with the conventional percolation, the 
difference is in the L dependence of H(p; L); H(p; L) is a 
decreasing function of L, therefore, there is no fixed point 
where the percolation probability IPs for all L intersect. 

As L becomes large, P5o(L) converges to a finite value 
p c and Ap(L) goes to zero as 

L -iM 

L -l/V2 



Pc 



P5o{L) 
Ap(L) 



(5) 
(6) 



with the exponents v\ and ~ 1/0.50, and the threshold 
p c sa 0.5221 (Fig. 3(c)). The corresponding value of /m, 
or f c , is /c/ffth ~ 0.926. The fact that v\ w v 2 suggests 
that there is a unique length scale £, which diverges as 



i ~ \P - Pc 



(7) 



with the exponent v « 2.0. Note that its value 2.0 is quite 
different from 4/3, or the value for the 2-d percolation. 9 

4-3 Order parameter 

The percolation order parameter P, or the ratio of the 
sites that belongs to the percolating cluster, is shown in 
Fig.4(a) for L = 64 - 8192. The value of P at p = p c is 
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Fig. 4. (a) The order parameter P vs the occupation ratio p for 
L = 64 ~ 8192. The vertical line is drawn at p = p c = 0.5221. 
(b)The scaling plot with (3 = 0.02 and u = 2.0. 

almost independent of L, thus the transition looks first 
order, but one can observe slight tendency of decrease in 
P at p = p c for the systems L > 1024, which implies the 
transition is of the second order with a small value of the 
exponent (3; 

P(p)~{p-Pcf- (8) 

The scaling plots 

P(p:L) = L-^f P ((p- Pc )L 1 /») (9) 

with the scaling function fp{x) are shown in Fig. 4(b) 
with (3 — 0.02 and v = 2.0. The plot looks reasonably 
good, but one may find some scatter. 

4-4 Cluster size distribution 

The cluster size distribution n(s) is introduced as the 
number of clusters with the size s divided by the number 
of lattice sites L 2 in the system, thus it is normalized as 

L 2 

J2sn{s)=p. (10) 

In Fig. 5, the cluster size distributions are shown in the 
logarithmic scale for p < p c (a) , p m p c (b) , and p > p c 
(c). In the case of p < p c (a), the distribution has finite 
width in the large L limit. As for the case of p > p c (c), 
it is also of finite width but with an isolated peak for 
each system size; The peak corresponds with the perco- 
lating clusters. The distribution at p ~ p c in Fig. 5(b) is 
intriguing; They show the power law distribution 

n{s)~ S - T (11) 

with the exponent r = 2.25, but there are peaks at the 
large cluster end of the distribution; The position of the 
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Fig. 5. The cluster size distributions below (a), at (b), and above 
(c) the percolation transition for various system sizes. The initial 
granular density is given by $m = 0.88(a), 0.926(b), and 0.95(c) 
with the resulting p ta 0.502(a), 0.522(b), and 0.531(c). The sys- 
tem sizes are L = 32 ~ 1024 (a), 128 ~ 4096 (b), and 128 ~ 1024 
(c). The other parameters are Af = 0.5/m, l/gth = 1 an d 
n s = 3. 

peak scales with almost L 2 which suggests they corre- 
sponds to the percolating cluster. This is consistent with 
the result of percolation probability in Fig. 3(b), which 
shows the system with any finite size is percolating at 

P = Pc- 

Fig. 6 shows the scaling plot 

n(s;L) = L- D \f n (s/L D ) (12) 

at p rts p c for D — 2 (a) and 1 (b) with r = 2.25 for the 
both. 

From Fig. 6 (a) with D = 2, one can see the peak posi- 
tion scales with the system area. This is consistent with 
the observation that the percolating clusters look com- 
pact and two-dimensional in Fig. 2. If one looks care- 
fully, however, one may find systematic deviation near 
the peak. This can be interpreted as the distribution de- 
viates from the power law (11) at s > s c (L), and s c (L) 
does not scale as L 2 but L D with another exponent D'. 

To see this, we plot the same data except for the 
two data points that correspond to the large percolat- 
ing clusters for each L; We tried various values of D', 
but D' = 1 seems to give least scatter around the master 




Fig. 6. The scaling plot of the cluster size distribution at the 
percolation transition p = 0.5221 (f M = 0.9259) for L = 128 ~ 
8192. (a) The data for s > 4 are plotted with D = 2 and r = 
2.25. (b) The data that correspond to the percolating clusters 
are removed in the plot with D' = 1 and r = 2.25, using the 
same scaling form (12). The solid lines indicate the line with the 
slope —2.25. 
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2-d percolation 
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Table I. The critical exponents at the percolation threshold for 
the sweeping interface model and the conventional 2-d percola- 
tion. 9 

curve(Fig.6(b)). The fact that D and D' are quite dif- 
ferent simply means that the boundary effect has strong 
influence on the cluster distribution, but the exponent 
D' should not be associated with the fractal dimension 
of the geometrical structure of clusters. 

5. Discussions and summary 

The exponents obtained above are listed with those 
for the two-dimensional conventional percolation transi- 
tion(Table I); They are clearly different, suggesting that 
the nature of the transition is different from that of the 
conventional percolation. This may not be very surpris- 
ing because the labyrinthine patterns analyzed here are 
results of the dynamical process of sweeping, which could 
induce long range correlation, while the conventional per- 
colation transition is for the pattern by random deposi- 
tion. 

Most remarkable difference, however, is the behavior 
of the percolation probability Ii(p;L) in Fig. 3(b). For 
the sweeping interface model, the H(p; L) is a decreasing 
function of L for any p, thus H(p; i)'s for different L do 
not intersect, while those for the percolation intersect at 
p = p c . This feature, however, may be explained by the 
fact that the clusters produced by the sweeping dynamics 
have branching structure with virtually no loops; Even 
in a very large cluster, the path that connects any given 
two sites in the cluster is essentially unique, 10 thus, for a 
fixed p, the percolating cluster is more difficult to appear 
in a larger system because breakage of a path at any 
site during drying process may break the cluster apart. 
One consequence of this feature is that the finite system 
is almost always percolated at the percolation threshold 
p = p c for the infinite system. 

Another aspect of the loopless structure of cluster is 
that the cluster structure of the system changes drasti- 
cally when the system is trimmed. Fig. 7 shows the cluster 
structure of the system with L = 1024 and its subsys- 
tems for p p c . The upper left plate is the original 
system of L — 1024. The largest percolating cluster col- 
ored by grey is dominating the system. The upper right, 
lower left, and lower right plates show the central parts of 
the original system of the size L = 512, 256, and 128, re- 
spectively. The cluster structure are analyzed only within 
the subsystems, i.e. the clusters that are connected only 
outside of the subsystems are disconnected in the subsys- 
tem. One can see there are no percolating clusters in the 
smallest subsystem because the original percolating clus- 
ter is fragmented into small pieces when we look at only 
a smaller part of it. This should be a general feature for 
the labyrinthine structure; There is virtually no loop in 
clusters, therefore, there is usually only one connecting 
path between any two sites in a cluster. Consequently, 
clusters are easily fragmented into pieces when they are 




Fig. 7. (Color on line) Cluster structure at the percolation thresh- 
old p = 0.522. The cluster structure of the whole system of the 
size L = 1024 (upper left), and its subsets at the central part of 
the size 512 (upper right), 256 (lower left), and 128 (lower right) 
are shown. 



trimmed to fit into a smaller part of the system. 

This may pose some problems in the analysis of ex- 
perimental data; In the experiment, it should be difficult 
to analyse whole systems because the regions near the 
boundary arc usually disturbed. In the finite size scaling 
analysis, the samples with smaller sizes are generated by 
trimming larger samples, 2 but the present analysis sug- 
gests that such a procedure changes the cluster structure 
and results may well depend upon the size of the original 
system. 

In summary, being stimulated by the labyrinthine pat- 
tern produced in the drying process, we have constructed 
the lattice model of the sweeping interface based on the 
invasion percolation, and demonstrated the model pro- 
duces similar patterns. The cluster analysis of the result- 
ing patterns shows that there is a percolation transition 
upon changing the density in the initial state, but the 
loopless structure of the labyrinthine patterns obtained 
by the sweeping dynamics make the nature of the tran- 
sition different from the conventional percolation. 
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